Dynamic Heterochromatin States in Anisotropic Nuclei of Cells on Aligned Nanofibers.

The cancer cell nucleus deforms as it invades the interstitial spaces in tissues and the tumor microenvironment. While alteration of the chromatin structure in a deformed nucleus is expected and documented, the chromatin structure in the nuclei of cells on aligned matrices has not been elucidated. In this work we elucidate the spatiotemporal organization of heterochromatin in the elongated nuclei of cells on aligned nanofibers with stimulated emission depletion nanoscopy and fluorescence correlation spectroscopy. We show that the anisotropy of nuclei is sufficient to drive H3K9me3-heterochromatin alterations, with enhanced H3K9me3 nanocluster compaction and aggregation states that otherwise are indistinguishable from diffraction-limited microscopy. We interrogated the higher-order heterochromatin structures within major chromatin compartments in anisotropic nuclei and discovered a wider spatial dispersion of nanodomain clusters in the nucleoplasm and condensed larger nanoclusters near the periphery and pericentromeric heterochromatin. Upon examining the spatiotemporal dynamics of heterochromatin in anisotropic nuclei, we observed reduced mobility of the constitutive heterochromatin mark H3K9me3 and the associated heterochromatin protein 1 (HP1α) at the nucleoplasm and periphery regions, correlating with increased viscosity and changes in gene expression. Since heterochromatin remodeling is crucial to genome integrity, our results reveal an unconventional H3K9me3 heterochromatin distribution, providing cues to an altered chromatin state due to perturbations of the nuclei in aligned fiber configurations.

Cell Fragment Formation, Migration, and Force Exertion on Extracellular Mimicking Fiber Nanonets.

Cell fragments devoid of the nucleus play an essential role in intercellular communication. Mostly studied on flat 2D substrates, their origins and behavior in native fibrous environments remain unknown. Here, cytoplasmic fragments' spontaneous formation and behavior in suspended extracellular matrices mimicking fiber architectures (parallel, crosshatch, and hexagonal) are described. After cleaving from the parent cell body, the fragments of diverse shapes on fibers migrate faster compared to 2D. Furthermore, while fragments in 2D are mostly circular, a higher number of rectangular and blob-like shapes are formed on fibers, and, interestingly, each shape is capable of forming protrusive structures. Absent in 2D, fibers' fragments display oscillatory migratory behavior with dramatic shape changes, sometimes remarkably sustained over long durations (>20 h). Immunostaining reveals paxillin distribution along fragment body-fiber length, while Forster Resonance Energy Transfer imaging of vinculin reveals mechanical loading of fragment adhesions comparable to whole cell adhesions. Using nanonet force microscopy, the forces exerted by fragments are estimated, and peculiarly small area fragments can exert forces similar to larger fragments in a Rho-associated kinase dependent manner. Overall, fragment dynamics on 2D substrates are insufficient to describe the mechanosensitivity of fragments to fibers, and the architecture of fiber networks can generate entirely new behaviors.

Author Correction: Force-exerting perpendicular lateral protrusions in fibroblastic cell contraction.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Force-exerting perpendicular lateral protrusions in fibroblastic cell contraction.

Aligned extracellular matrix fibers enable fibroblasts to undergo myofibroblastic activation and achieve elongated shapes. Activated fibroblasts are able to contract, perpetuating the alignment of these fibers. This poorly understood feedback process is critical in chronic fibrosis conditions, including cancer. Here, using fiber networks that serve as force sensors, we identify "3D perpendicular lateral protrusions" (3D-PLPs) that evolve from lateral cell extensions named twines. Twines originate from stratification of cyclic-actin waves traversing the cell and swing freely in 3D to engage neighboring fibers. Once engaged, a lamellum forms and extends multiple secondary twines, which fill in to form a sheet-like PLP, in a force-entailing process that transitions focal adhesions to activated (i.e., pathological) 3D-adhesions. The specific morphology of PLPs enables cells to increase contractility and force on parallel fibers. Controlling geometry of extracellular networks confirms that anisotropic fibrous environments support 3D-PLP formation and function, suggesting an explanation for cancer-associated desmoplastic expansion.

ECM in Differentiation: A Review of Matrix Structure, Composition and Mechanical Properties.

Stem cell regenerative potential owing to the capacity to self-renew as well as differentiate into other cell types is a promising avenue in regenerative medicine. Stem cell niche not only provides physical scaffolding but also possess instructional capacity as it provides a milieu of biophysical and biochemical cues. Extracellular matrix (ECM) has been identified as a major dictator of stem cell lineage, thus understanding the structure of in vivo ECM pertaining to specific tissue differentiation will aid in devising in vitro strategies to improve the differentiation efficiency. In this review, we summarize details about the native architecture, composition and mechanical properties of in vivo ECM of the early embryonic stages and the later adult stages. Native ECM from adult tissues categorized on their origin from respective germ layers are discussed while engineering techniques employed to facilitate differentiation of stem cells into particular lineages are noted. Overall, we emphasize that in vitro strategies need to integrate tissue specific ECM biophysical cues for developing accurate artificial environments for optimizing stem cell differentiation.

Bioenergetics underlying single-cell migration on aligned nanofiber scaffolds.

Cell migration is centrally involved in a myriad of physiological processes, including morphogenesis, wound healing, tissue repair, and metastatic growth. The bioenergetics that underlie migratory behavior are not fully understood, in part because of variations in cell culture media and utilization of experimental cell culture systems that do not model physiological connective extracellular fibrous networks. In this study, we evaluated the bioenergetics of C2C12 myoblast migration and force production on fibronectin-coated nanofiber scaffolds of controlled diameter and alignment, fabricated using a nonelectrospinning spinneret-based tunable engineered parameters (STEP) platform. The contribution of various metabolic pathways to cellular migration was determined using inhibitors of cellular respiration, ATP synthesis, glycolysis, or glucose uptake. Despite immediate effects on oxygen consumption, mitochondrial inhibition only modestly reduced cell migration velocity, whereas inhibitors of glycolysis and cellular glucose uptake led to striking decreases in migration. The migratory metabolic sensitivity was modifiable based on the substrates present in cell culture media. Cells cultured in galactose (instead of glucose) showed substantial migratory sensitivity to mitochondrial inhibition. We used nanonet force microscopy to determine the bioenergetic factors responsible for single-cell force production and observed that neither mitochondrial nor glycolytic inhibition altered single-cell force production. These data suggest that myoblast migration is heavily reliant on glycolysis in cells grown in conventional media. These studies have wide-ranging implications for the causes, consequences, and putative therapeutic treatments aimed at cellular migration.

Inositol polyphosphate multikinase is a metformin target that regulates cell migration.

Metformin has been shown to alter cell adhesion protein expression, which is thought to play a role in its observed antitumor properties. We found that metformin treatment down-regulated integrin ß1 concomitant with the loss of inositol polyphosphate multikinase (IPMK) in murine myocytes, adipocytes, and hepatocytes. To determine if IPMK was upstream of integrin ß1 expression, we examined IPMK-/- mouse embryonic fibroblast cells and found that integrins ß1 and ß3 gene expression was reduced by half, relative to wild-type cells, whereas focal adhesion kinase (FAK) activity and Rho/Rac/Cdc42 protein levels were increased, resulting in migration defects. Using nanonet force microscopy, we determined that cell:extracellular matrix adhesion and cell contractility forces were decreased, confirming the functional relevance of integrin and Rho protein dysregulation. Pharmacological studies showed that inhibition of both FAK1 and proline-rich tyrosine kinase 2 partially restored integrin ß1 expression, suggesting negative regulation of integrin ß1 by FAK. Together our data indicate that IPMK participates in the regulation of cell migration and provides a potential link between metformin and wound healing impairment.-Tu-Sekine, B., Padhi, A., Jin, S., Kalyan, S., Singh, K., Apperson, M., Kapania, R., Hur, S. C., Nain, A., Kim, S. F. Inositol polyphosphate multikinase is a metformin target that regulates cell migration.

Construction of PRDM9 allele-specific recombination maps in cattle using large-scale pedigree analysis and genome-wide single sperm genomics.

PRDM9 contributes to hybrid sterility and species evolution. However, its role is to be confirmed in cattle, a major domesticated livestock species. We previously found an association near PRDM9 with cattle recombination features, but the causative variants are still unknown. Using millions of genotyped cattle with pedigree information, we characterized five PRDM9 alleles and generated allele-specific recombination maps. By examining allele-specific recombination patterns, we observed the impact of PRDM9 on global distribution of recombination, especially in the two ends of chromosomes. We also showed strong associations between recombination hotspot regions and functional mutations within PRDM9 zinc finger domain. More importantly, we found one allele of PRDM9 to be very different from others in both protein composition and recombination landscape, indicating the causative role of this allele on the association between PRDM9 and cattle recombination. When comparing recombination maps from sperm and pedigree data, we observed similar genome-wide recombination patterns, validating the quality of pedigree-based results. Collectively, these evidence supported PRDM9 alleles as causal variants for the reported association with cattle recombination. Our study comprehensively surveyed the bovine PRDM9 alleles, generated allele-specific recombination maps, and expanded our understanding of the role of PRDM9 on genome distribution of recombination.

Aligned fibers direct collective cell migration to engineer closing and nonclosing wound gaps.

Cell emergence onto damaged or organized fibrous extracellular matrix (ECM) is a crucial precursor to collective cell migration in wound closure and cancer metastasis, respectively. However, there is a fundamental gap in our quantitative understanding of the role of local ECM size and arrangement in cell emergence-based migration and local gap closure. Here, using ECM-mimicking nanofibers bridging cell monolayers, we describe a method to recapitulate and quantitatively describe these in vivo behaviors over multispatial (single cell to cell sheets) and temporal (minutes to weeks) scales. On fiber arrays with large interfiber spacing, cells emerge (invade) either singularly by breaking cell-cell junctions analogous to release of a stretched rubber band (recoil), or in groups of few cells (chains), whereas on closely spaced fibers, multiple chains emerge collectively. Advancing cells on fibers form cell streams, which support suspended cell sheets (SCS) of various sizes and curvatures. SCS converge to form local gaps that close based on both the gap size and shape. We document that cell stream spacing of 375 µm and larger hinders SCS advancement, thus providing abilities to engineer closing and nonclosing gaps. Altogether we highlight the importance of studying cell-fiber interactions and matrix structural remodeling in fundamental and translational cell biology.

Characterization of genome-wide segmental duplications reveals a common genomic feature of association with immunity among domestic animals.

BACKGROUND: Segmental duplications (SDs) commonly exist in plant and animal genomes, playing crucial roles in genomic rearrangement, gene innovation and the formation of copy number variants. However, they have received little attention in most livestock species. RESULTS: Aiming at characterizing SDs across the genomes of diverse livestock species, we mapped genome-wide SDs of horse, rabbit, goat, sheep and chicken, and also enhanced the existing SD maps of cattle and pig genomes based on the most updated genome assemblies. We adopted two different detection strategies, whole genome analysis comparison and whole genome shotgun sequence detection, to pursue more convincing findings. Accordingly we identified SDs for each species with the length of from 21.7 Mb to 164.1 Mb, and 807 to 4,560 genes were harboured within the SD regions across different species. More interestingly, many of these SD-related genes were involved in the process of immunity and response to external stimuli. We also found the existence of 59 common genes within SD regions in all studied species except goat. These common genes mainly consisted of both UDP glucuronosyltransferase and Interferon alpha families, implying the connection between SDs and the evolution of these gene families. CONCLUSIONS: Our findings provide insights into livestock genome evolution and offer rich genomic sources for livestock genomic research.

Ruminant-specific multiple duplication events of PRDM9 before speciation.

BACKGROUND: Understanding the genetic and evolutionary mechanisms of speciation genes in sexually reproducing organisms would provide important insights into mammalian reproduction and fitness. PRDM9, a widely known speciation gene, has recently gained attention for its important role in meiotic recombination and hybrid incompatibility. Despite the fact that PRDM9 is a key regulator of recombination and plays a dominant role in hybrid incompatibility, little is known about the underlying genetic and evolutionary mechanisms that generated multiple copies of PRDM9 in many metazoan lineages. RESULTS: The present study reports (1) evidence of ruminant-specific multiple gene duplication events, which likely have had occurred after the ancestral ruminant population diverged from its most recent common ancestor and before the ruminant speciation events, (2) presence of three copies of PRDM9, one copy (lineages I) in chromosome 1 (chr1) and two copies (lineages II & III) in chromosome X (chrX), thus indicating the possibility of ancient inter- and intra-chromosomal unequal crossing over and gene conversion events, (3) while lineages I and II are characterized by the presence of variable tandemly repeated C2H2 zinc finger (ZF) arrays, lineage III lost these arrays, and (4) C2H2 ZFs of lineages I and II, particularly the amino acid residues located at positions -1, 3, and 6 have evolved under strong positive selection. CONCLUSIONS: Our results demonstrated two gene duplication events of PRDM9 in ruminants: an inter-chromosomal duplication that occurred between chr1 and chrX, and an intra-chromosomal X-linked duplication, which resulted in two additional copies of PRDM9 in ruminants. The observation of such duplication between chrX and chr1 is rare and may possibly have happened due to unequal crossing-over millions of years ago when sex chromosomes were independently derived from a pair of ancestral autosomes. Two copies (lineages I & II) are characterized by the presence of variable sized tandem-repeated C2H2 ZFs and evolved under strong positive selection and concerted evolution, supporting the notion of well-established Red Queen hypothesis. Collectively, gene duplication, concerted evolution, and positive selection are the likely driving forces for the expansion of ruminant PRDM9 sub-family.

Positive Selection Drives Rapid Evolution of the meq Oncogene of Marek's Disease Virus.

Marek's disease (MD), caused by Marek's disease virus (MDV), a poultry-borne alphaherpesvirus, is a devastating disease of poultry causing an estimated annual loss of one billion dollars to poultry producers, worldwide. Despite decades of control through vaccination, MDV field strains continue to emerge having increased virulence. The evolutionary mechanism driving the emergence of this continuum of strains to increased MDV virulence, however, remains largely enigmatic. Increase in MDV virulence has been associated with specific amino acid changes within the C-terminus domain of Mareks's EcoRI-Q (meq)-encoded oncoprotein. In this study, we sought to determine whether the meq gene has evolved adaptively and whether past vaccination efforts have had any significant effect on the reduction or increase of MDV diversity over time. Our analysis suggests that meq is estimated to be evolving at a much faster rate than most dsDNA viruses, and is comparable with the evolutionary rate of RNA viruses. Interestingly, most of the polymorphisms in meq gene appear to have evolved under positive selection and the time of divergence at the meq locus coincides with the period during which the poultry industry had undergone transitions in management practices including the introduction and widespread use of live attenuated vaccines. Our study has revealed that the decades-long use of vaccines did not reduce MDV diversity, but rather had a stimulating effect on the emergence of field strains with increased genetic diversity until the early 2000s. During the years 2004-2005, there was an abrupt decline in the genetic diversity of field isolates followed by a recovery from this bottleneck in the year 2010. Collectively, these data suggest that vaccination seems to not have had any effect on MDV eradication, but rather had a stimulating effect on MDV emergence through adaptation.

Evidence for Within-Host Genetic Recombination among the Human Pegiviral Strains in HIV Infected Subjects.

The non-pathogenic Human Pegivirus (HPgV, formerly GBV-C/HGV), the most prevalent RNA virus worldwide, is known to be associated with reduced morbidity and mortality in HIV-infected individuals. Although previous studies documented its ubiquity and important role in HIV-infected individuals, little is known about the underlying genetic mechanisms that maintain high genetic diversity of HPgV within the HIV-infected individuals. To assess the within-host genetic diversity of HPgV and forces that maintain such diversity within the co-infected hosts, we performed phylogenetic analyses taking into account 229 HPgV partial E1-E2 clonal sequences representing 15 male and 8 female co-infected HIV patients from Hubei province of central China. Our results revealed the presence of eleven strongly supported clades. While nine clades belonged to genotype 3, two clades belonged to genotype 2. Additionally, four clades that belonged to genotype 3 exhibited inter-clade recombination events. The presence of clonal sequences representing multiple clades within the HIV-infected individual provided the evidence of co-circulation of HPgV strains across the region. Of the 23 patients, six patients (i.e., five males and one female) were detected to have HPgV recombinant sequences. Our results also revealed that while male patients shared the viral strains with other patients, viral strains from the female patients had restricted dispersal. Taken together, the present study revealed that multiple infections with divergent HPgV viral strains may have caused within-host genetic recombination, predominantly in male patients, and therefore, could be the major driver in shaping genetic diversity of HPgV.

Life History of the Oldest Lentivirus: Characterization of ELVgv Integrations in the Dermopteran Genome.

Endogenous retroviruses are genomic elements formed by germline infiltration by originally exogenous viruses. These molecular fossils provide valuable information about the evolution of the retroviral family. Lentiviruses are an extensively studied genus of retroviruses infecting a broad range of mammals. Despite a wealth of information on their modern evolution, little is known about their origins. This is partially due to the scarcity of their endogenous forms. Recently, an endogenous lentivirus, ELVgv, was discovered in the genome of the Malayan colugo (order Dermoptera). This represents the oldest lentiviral evidence available and promises to lead to further insights into the history of this genus. In this study, we analyzed ELVgv integrations at several genomic locations in four distinct colugo specimens covering all the extant dermopteran species. We confirmed ELVgv integrations in all the specimens examined, which implies that the virus originated before the dermopteran diversification. Using a locus-specific dermopteran substitution rate, we estimated that the proviral integrations occurred 21-40 Ma. Using phylogenetic analysis, we estimated that ELVgv invaded an ancestor of today's Dermoptera in an even more distant past. We also provide evidence of selective pressure on the TRIM5 antiviral restriction factor, something usually taken as indirect evidence of past retroviral infections. Interestingly, we show that TRIM5 was under strong positive selection pressure only in the common dermopteran ancestor, where the ELVgv endogenization occurred. Further experiments are required to determine whether ELVgv participated in the TRIM5 selection.

Testis-enriched heat shock protein A2 (HSPA2): Adaptive advantages of the birds with internal testes over the mammals with testicular descent.

The molecular chaperone heat shock protein A2 (HSPA2), a member of the 70 kDa heat shock protein (HSP70) family, plays an important role in spermatogenesis and male fertility. Although HSPA2 is evolutionarily highly conserved across the metazoan lineages, the observation of striking differences in temperature-sensitive expressions, testicular physiology, spermatogenesis, as well as its role in male fertility indicates that avian and mammalian HSPA2 may exhibit distinct evolutionary trajectory. The present study reports that while mammalian HSPA2 is constrained by intense purifying selection, avian HSPA2 has been subjected to positive selection. The majority of the positively selected amino acid residues fall on the α-helix and ß-sheets of the peptide-binding domain located at the carboxyl-terminal region of the avian HSPA2. The detection of positively selected sites at the helix and ß-sheets, which are less tolerant to molecular adaptation, indicates an important functional consequence and contribution to the structural and functional diversification of the avian HSPA2. Collectively, avian HSPA2 may have an adaptive advantage over the mammals in response to heat stress, and therefore, mammals with testicular descent may be at a greater risk in the event of scrotal temperature rise.

Cattle Sex-Specific Recombination and Genetic Control from a Large Pedigree Analysis.

Meiotic recombination is an essential biological process that generates genetic diversity and ensures proper segregation of chromosomes during meiosis. From a large USDA dairy cattle pedigree with over half a million genotyped animals, we extracted 186,927 three-generation families, identified over 8.5 million maternal and paternal recombination events, and constructed sex-specific recombination maps for 59,309 autosomal SNPs. The recombination map spans for 25.5 Morgans in males and 23.2 Morgans in females, for a total studied region of 2,516 Mb (986 kb/cM in males and 1,085 kb/cM in females). The male map is 10% longer than the female map and the sex difference is most pronounced in the subtelomeric regions. We identified 1,792 male and 1,885 female putative recombination hotspots, with 720 hotspots shared between sexes. These hotspots encompass 3% of the genome but account for 25% of the genome-wide recombination events in both sexes. During the past forty years, males showed a decreasing trend in recombination rate that coincided with the artificial selection for milk production. Sex-specific GWAS analyses identified PRDM9 and CPLX1 to have significant effects on genome-wide recombination rate in both sexes. Two novel loci, NEK9 and REC114, were associated with recombination rate in both sexes, whereas three loci, MSH4, SMC3 and CEP55, affected recombination rate in females only. Among the multiple PRDM9 paralogues on the bovine genome, our GWAS of recombination hotspot usage together with linkage analysis identified the PRDM9 paralogue on chromosome 1 to be associated in the U.S. Holstein data. Given the largest sample size ever reported for such studies, our results reveal new insights into the understanding of cattle and mammalian recombination.

A testis-specific gene within a widely expressed gene: Contrasting evolutionary patterns of two differentially expressed mammalian proteins encoded by a single gene, CAMK4.

Understanding the patterns of genetic variations within fertility-related genes and the evolutionary forces that shape such variations is crucial in predicting the fitness landscapes of subsequent generations. This study reports distinct evolutionary features of two differentially expressed mammalian proteins [CaMKIV (Ca(2+) /calmodulin-dependent protein kinase IV) and CaS (calspermin)] that are encoded by a single gene, CAMK4. The multifunctional CaMKIV, which is expressed in multiple tissues including testis and ovary, is evolving at a relatively low rate (0.46-0.64 × 10(-9) nucleotide substitutions/site/year), whereas the testis-specific CaS gene, which is predominantly expressed in post-meiotic cells, evolves at least three to four times faster (1.48-1.98 × 10(-9) substitutions/site/year). Concomitantly, maximum-likelihood-based selection analyses revealed that the ubiquitously expressed CaMKIV is constrained by intense purifying selection and, therefore, remained functionally highly conserved throughout the mammalian evolution, whereas the testis-specific CaS gene is under strong positive selection. The substitution rates of different mammalian lineages within both genes are positively correlated with GC content, indicating the possible influence of GC-biased gene conversion on the estimated substitution rates. The observation of such unusually high GC content of the CaS gene (≈74%), particularly in the lineage that comprises the bovine species, suggests the possible role of GC-biased gene conversion in the evolution of CaS that mimics positive selection.

Time-dependent selection pressure on two arthropod-borne RNA viruses in the same serogroup.

Understanding the genetic basis of viral adaptation to taxonomically diverse groups of host species inhabiting different eco-climatic zones is crucial for the discovery of factors underpinning the successful establishment of these infectious pathogens in new hosts/environments. To gain insights into the dynamics of nonsynonymous (dN) and synonymous substitutions (dS) and the ratio between the two (ω=dN/dS), we analyzed the complete nucleotide coding sequence data of the M segment, which encodes glycoproteins of two negative-sense RNA viruses, Akabane virus (AKV) and Schmallenberg virus (SBV) that belong to the same serogroup. While AKV is relatively older and has been circulating in ruminant populations since 1970s, SBV was first reported in 2011. The ω was estimated to be 1.67 and 0.09 for SBV and AKV, respectively, and the estimated mutation rate of SBV is at least 25 times higher than that of AKV. Given the different evolutionary stages of the two viruses, most of the slightly deleterious mutations were likely purged out or kept in low frequency in the AKV genome, whereas positive selection together with the accumulation of slightly deleterious mutations might contribute to such an inflated mutation rate of SBV. The evolutionary distance (d) is nonlinearly and negatively correlated with ω, but is positively correlated with dN and dS. Collectively, the different patterns in ω, dN, dS, and d between AKV and SBV identified in this study provide empirical evidence for a time-dependent selection pressure.

Genetic and epidemiological insights into the emergence of peste des petits ruminants virus (PPRV) across Asia and Africa.

Small ruminants are important components in the livelihood of millions of households in many parts of the world. The spread of the highly contagious peste des petits ruminants (PPR) disease, which is caused by an RNA virus, PPRV, across Asia and Africa remains a major concern. The present study explored the evolutionary and epidemiological dynamics of PPRV through the analyses of partial N-gene and F-gene sequences of the virus. All the four previously described PPRV lineages (I-IV) diverged from their common ancestor during the late-19(th) to early-20(th) century. Among the four lineages, PPRV-IV showed pronounced genetic structuring across the region; however, haplotype sharing among the geographic regions, together with the presence of multiple genetic clusters within a country, indicates the possibility of frequent mobility of the diseased individuals across the region. The gradual decline in the effective number of infections suggests a limited genetic variation, which could be attributed to the effective vaccination that has been practiced since 1990s. However, the movement of infected animals across the region likely contributes to the spread of PPRV-IV. No evidence of positive selection was identified from this study.

Molecular evolutionary and epidemiological dynamics of genotypes 1G and 2B of rubella virus.

Rubella Virus (RV), which causes measles-like rashes in children, puts millions of infants at risk of congenital defects across the globe. Employing phylogenetic approaches to the whole genome sequence data and E1 glycoprotein sequence data, the present study reports the substitution rates and dates of emergence of all thirteen previously described rubella genotypes, and gains important insights into the epidemiological dynamics of two geographically widely distributed genotypes 1G and 2B. The overall nucleotide substitution rate of this non-vector-borne RV is in the order of 10-3 substitutions/site/year, which is considerably higher than the substitution rates previously reported for the vector-borne alphaviruses within the same family. Currently circulating strains of RV share a common ancestor that existed within the last 150 years, with 95% Highest Posterior Density values ranging from 1868 to 1926 AD. Viral strains within the respective genotypes began diverging between the year 1930 s and 1980 s. Both genotype 1G and 2B have shown a decline in effective number of infections since 1990 s, a period during which mass immunization programs against RV were adapted across the globe. Although both genotypes showed some extent of spatial genetic structuring, the analyses also depicted an inter-continental viral dispersal. Such a viral dispersal pattern could be related to the migration of infected individuals across the regions coupled with a low coverage of MMR vaccination.